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TURBO DECODING WITH ITERATIVE ESTIMATION OF CHANNEL PARAMETERS 

Field of the invention 

The present invention relates to a method for decoding a convolutionally coded 
5 input data signal y comprising multiplying the input data signal with a scaling factor 
, demultiplexing the multiplied input data signal L^y , e.g. in parity signals and a 
systematic signal, and turbo decoding the demultiplexed input data signal L^S in order 

to obtain decoder output likelihood ratio data. In a further aspect, the present invention 
relates to a decoder device for decoding a convolutionally coded input data signal y 
10 comprising a multiplication element for multiplying a received input data signal y with 
a scaling factor , a demultiplexer for demultiplexing the multiplied input data signal 
L^y , e.g. in parity signals and a systematic signal, and a turbo decoder for decoding the 
demultiplexed input data signal L^S in order to obtain decoder output likelihood ratio 
data. 

15 

Prior art 

Such a method of decoding data and a decoder device are known from American 
patent publication US-B-6,393,0765 which describes a method for decoding turbo codes 
using data scaling. Convolutionally coded input data is decoded in a near ideal manner. 
20 A portion of the input data is buffered, after which a mean of the data in the portion of 
the input data is calculated. Then, a root-mean-square value of the portion is calculated 
using the mean. A scaling factor is derived from the root-mean-square value, which 
scaling factor is used to scale the portion of the input data before the turbo decoding 
step. 

25 This disclosed method has the disadvantage that the scaling factor computation is 

not based on a-posteriori likelihood information, which leads to additional loss. 
Especially in mobile applications using these kind of coding, every additional loss has a 
negative effect on system performance. 

The publication of El-Gamal, H., High capacity Synclironous FH/SSMA 

30 networks with turbo coding, IEEE Intern. Symp. on Spread Spectrum Tecliniques, 
1998, p973-977, provides a similar method for scaling the input data. The method 
consists of computing the noise variation, which is subsequently applied to form log- 
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likelihood values. This method has the disadvantage that it assxunes the signal 
amplitude as nomialized to unity. This is a simplification not valid for practical 
receiving equipment, and therefore negatively impacts the system performance. 

Summary of the invention 

The present invention seeks to provide an improved decoding scheme for use in 
transceivers using a maximum-a-posteriori (MAP) decoder, or related techniques such 
as Logarithmic MAP (LOGMAP). 

According to the present invention, a method according to the preamble defined 
above is provided, in which the scaling factor is updated for a next iteration in 
dependence on a combination of a posteriori likelihood data based on turbo decoded 
output data A and a priori likelihood data based on the demultiplexed signal L^S . By 
combining these a posteriori and a priori likelihood data, it is possible to improve the 
perforaiance of the decoding method without the need of much additional hardware or 
software resources. 

In an embodiment of the present invention the scaling factor is updated using an 
estimate of the mean value of the signal amplitude c and an estimate of the noise 
variation . The update of the scaling factor can e.g. be calculated according to 

L = — vL^ , in which is the updated scaling factor. The estimate of the mean 

value of the signal amplitude may be equal to c = — ^ sgn(A, ) ^L^-s^ , 



where N is the number of bits in a coding block of the input data signal, is the f' 
systematic bit, c is the estimation of the amplitude of the scaled systematic bits • 
and A,, is the log-likelihood ratio resulting fi:om the most recent turbo decoder 
iteration. 

In a further embodiment, the noise variance estimation may equal 




1 



2 (si - 1) ^ . (1) + (si ^if^PM-K; 



the probability of the 



z^'bit being zero is estimated like Pr{:x:,. = o} = (0) = 




wo 2005/062473 



PCT/NL2003/000924 



3 

and ttie probability of that bit being one like Pr{x^ = l} = (1) = — = 1 - (0) ; 

l + e ' 

L 'S 

Ihe normalised systematic bits s'^ are calculated as = ' ; 

c 

and where is a bias correction computed as K = — 5] ^ " (^) " ^ 0)) 

5 This makes the present method accurate and fast, so that it secures the advantage 

(simulations show the loss compared to the theoretical optimum is about 0.03dB), and 
is applicable in fast fading or time-variant channels. It has been shown that about 0.1 to 
0.2dB in sensitivity improvement for turbo codes can be achieved, compared to prior 
art methods. 

10 In a further embodiment, the scaling factor is initialized either as a fixed 

value, the result of an initial number of iterations using a known algorithm, filtering 
over subsequent iterations and coding blocks, or SNR/SIR estimation at the mput data 
signal. Initialization can thus be accomplished using very simple solutions or more 
complex but well known solutions. 

15 In an even further embodiment of the present invention, the method further 

comprises calculating the variation of the scaling factor in subsequent iterations and, 
when the variation after a predetermined number of iterations is above a predetermined 
threshold value, reverting to a different scaling factor calculation method and/or turbo 
decoding method. The different scaling factor calculation method may e.g. be a fixed 

20 scaling factor, a Log-MAX method or a SOVA method (Soft Output Viterbi 

Algorithm). In this manner, a simple divergence monitoring of the iterative method can 
be accomplished, using known algorithms as back-up. 

In a further aspect of the present invention, a decoder device is provided 
according to the preamble defined above, in which the decoder device further 

25 comprises an adaptive scaling element which is arranged to update the scaling factor 
for a next iteration based on a combination of a posteriori likelihood data based on 
turbo decoded output data and a priori likelihood data based on the demultiplexed 
signal. 

The adaptive scaling element may be further arranged to execute the present 
30 method. Accordingly, the present decoder device provides an improved performance 
over prior art devices. Because the sensitivity improvement only requires a small 
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hardware cost, the sensitivity improvement is a pure bonus. The important system 
aspects, which influence sensitivity, are signal strength, i.e. range/coverage/battery life, 
noise figure, and interference, i.e., air interface capacity. Hence, the advantage of the 
present invention can be translated mto a coverage increase (<1% more range), 
5 mcreased mobile battery life time (2 to 4% less transmit power), relaxed noise figure 
requirements (0.2dB less), or improved capacity users on the air interface (2 to 4% 
more). This yields a larger cell-coverage and lower power usage for the mobile, which 
in turn gives less interference, and therefore more capacity. 

In an even further aspect, the present invention provides a computer program 
10 product, which comprises computer executable code, which when loaded on a 

processing system, provides the processmg system with the capability to execute the 
present method. The processing system may comprise a microprocessor and peripheral 
equipment, a digital signal processor or a combination of both to execute the present 
method. 



15 



Short description of drawings 

The present invention will be discussed in more detail below, using a nximber of 
exemplary embodiments, with reference to the attached drawings, in which 

Fig. 1 shows a block diagram of an embodiment of the decoding scheme 
20 according to the present invention. 



Detailed description of exemplary embodiments 

The present application may be advantageously applied in decoding of systematic 
forward error-correcting codes, such as the parallel concatenated convolutional turbo 
25 codes found in third generation mobile telecommunication systems (3 GPP). The 
present invention can be employed in transceivers using a maximum-a-posteriori 
(MAP) decoder, or related techniques, such as Log-MAP. 

Iterative MAP and Log-MAP decoding are asympotically optimum. Sub-optimal 
decoding algorithms, such as soft output Viterbi algorithm (SOVA) or approximated 
30 LOGMAP (Log-Max), are simpler implementation-wise, yet about 0.5dB short of the 
MAP and Log-MAP performance. 

The problem with achieving the full extent of MAP performance is the 
requirement that the input data is defmed as a (log) likelihood ratio. With the present 



wo 2005/062473 



PCT/NL2003/000924 



invention it is possible to scale the decoder input data in a time-variant, fading channel, 
quickly and accurately, in order to form (log)likelihood values. 

The present method uses the a-posteriori likelihood data to form probability 
information of the decoded bits. This is combined with the decoder a-priori likelihood 
5 data, in order to optimally determine the desired signal mean and noise variance. These 
are used to iteratively scale the a-priori information. 

The scaling factor should be initialized. The invention foresees in several 
different methods for this: 

Use the scaling factor from the previous coding block (filtered or not); 
10 - Run the first iteration(s) with SOVA or Log-Max; 

Deploy any method as initial value setting; or 

Any combination of these. 

Fig. 1 shows a block diagram of the iterative scaling algorithm 10 for decoding 
convolutionally coded data. The received symbols y are multiplied by a scaling factor 

15 using a multiplier 8 before they are de-multiplexed using a demultiplexer 6 into 

systematic bits s. and parity bits par^ , par^ that are input to a turbo decoder 5. After 
an iteration of the decoder 5 an adaptive scaler 7 makes an estimation of the signal 
amplitude and the noise variance, based on the decoder output log-likelihood ratios A^. . 
From these estimations a new scaling factor is computed for the next iteration. 

20 The turbo decoder 5 comprises as shown in Fig. 1, interleavers 13, 14, de- 

interleaver 15, and decoders 1 1, 12. More in detail, the turbo decoder 5 comprises jfirst 
and second soft input soft output (SISO) decoder sections 1 1, 12 which receive data 
from the demultiplexer 6. The fimctioning of the turbo decoder 5 is known to the 
person skilled in the art (see e.g. US-B-6,393,076) and needs no further explanation in 

25 this description. 

In the following the operation of the adaptive scaler 7 is explained in detail. The 
calculations performed by the adaptive scaler 7 are described as algorithms, 
calculations and the like. 

For the person skilled in the art, it will be clear that the algorithms and 

30 calculations as performed by the various blocks and elements of the decoder algorithm 
10 may be implemented using software, hardware resoxxrces or a combination of both, 
such as analogue electronics, logical circuitry, signal processors, etc. 



wo 2005/062473 PCT/NL2003/000924 

6 

The signal amplitude is computed in block 20 from the log-likelihood ratios A,, 
resulting from the most recent turbo decoder 5 iteration, and from the systematic bits as 
follows: 

5 

where N is the number of bits in a coduig block and is the systematic bit. c is 
the estimation of the amplitude of the scaled systematic bits L^-s^; if is the 

estimated symbol energy of the signal component of input sequence , then ^JM^ is 
the estimated amplitude of this signal component; then we fibnd that 

10 c^L^'4^s => ^^=[t] 

In block 21 the estimated ampUtude c is used to normalize the systematic bits s'^ 



givmg: 



From the log-likelihood ratios A,, an estimation of the logical bit probabilities 
1 5 can be computed in block 22; the probability of the bit being 0 is esthnated like 

Pr{x,=0} = P,(0)=^-^ 

and the probabiUty of that bit being 1 like 

Pr{x, = 1} = (1) = = l-P, (0) 

20 

The non-linear relationship between A,, and logical bit probabilities Pi(0) and 
Pi(l) can be computed or calculated based on a look-up-table. 

These probabilities together with the normalised systematic bits s'^ can be used to 

compute an estimation o*,^, of the variance of the noise in block 23, according to 
25 ■?,(!) + (si +ir'PM-K 

where K is sl bias correction computed as follows: 
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Then the relation between cj^, and and is given by 

When the above equations are rewritten to 
and use 

c 



the new estimated optimum scaling factor can be computed in block 24 as: 

1=4.- 



N, 2E,ai, c-ai " 

10 Summarizing, the new scaling factor in the amplitude and variance estimations, and in 
the current scaling factor can be expressed as follows: 



This calculation is performed using multiplier 25 and memory latch 26, after 
15 which the current scaling factor is applied again to the input datay. 

Since the turbo decoder 5 is iterative itself, merging it with another iterative 
algorithm (that of the adaptive scaler 7) requires careful consideration in order to avoid 
divergence or oscillation. Initialisation of the scaling factor to a reasable starting 

value proves to be sttfficient In the absence of proper initialisation information, the 
20 first decoding iterations can be done with known LogMax or SOVA techniques. These 
are not scaling factor dependent, and can therefore be used to let the a-posteriori 
likelihood values converge. Other techniques to acquire initialisation information are 
filtering over subsequent iterations and coding blocks, and SNR/SIR estimation at the 
input, see for instance the method described in Sxxmmers, T.A., Wilson, S.G., SNR 
25 Mismatch and Online Estimation in Turbo Decoding, IEEE Tr. On COM, Vol. 46, 
No.4., April 1998, pp. 421-423. 

Moreover, divergence can be detected by monitoring the development of the 
scaling factor . Normally, it should converge in a few iterations, and then hardly 
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change. In the case of detected divergence, one can fix the scaling factor , or switch 

to Log-Max or SOVA. 

The description of the embodiment above, and set forth in the claims, is based on 
the LOG-MAP type of turbo decoder. A person skilled in the art will recognize that the 
5 ideas of the disclosed method are not restricted to LOG-MAP decoding, but can also be 
formulated for variations on the LOG-MAP decoding method, and for MAP decoding. 



